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Abstract. Given a set of wireless links, a fundamental problem is to find the largest subset that 
can transmit simultaneously, within the SINR model of interference. Significant progress on this 
problem has been made in recent years. In this note, we study the problem in the setting where 
we are given a fixed set of arbitrary powers each sender must use, and an arbitrary gain matrix 
defining how signals fade. This variation of the problem appears immune to most algorithmic 
approaches studied in the literature. Indeed it is very hard to approximate since it generalizes the 
max independent set problem. Here, we propose a simple semi-definite programming approach to 
the problem that yields constant factor approximation, if the optimal solution is strictly larger than 



half of the input size. 
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1. Introduction 



Q . We consider the fundamental problem of wireless network capacity. Given is a set L = {l\,l2, ■ ■ ■ An} 

of links, where each link £ v represents a communication request from a sender s v to a receiver r v . 
We are also given, for every £ v £ L, a transmission power P v > 0. The powers received from senders 
to receivers are defined by an n x n dimensional gain matrix G with positive entries. Specifically, 
the signal received from s v at r w is G wv ■ P v . Thus an instance in this model can be described by 
the tuple (L, P, G) where P is the vector of the power assignments P v for all £ v . 
1/-) Simultaneously communicating links interfere with each other, following the physical model or 

"SINR model" of interference. Due to its higher fidelity to reality [101 [22| 126] . this model of 
interference has recently gained substantial attention in the analysis of wireless networks. In this 
model, a receiver r v successfully receives a message from a sender s v if and only if the following 
condition holds: 

i-\ \ G vv ■ P v 

(!) v r p~^n - /3 ' 

where N is a universal constant denoting the ambient noise, f3 > 1 denotes the minimum SINR 
(signal-to-interference-noise-ratio) required for a message to be successfully received, and S is the 
set of concurrently scheduled links in the same slot. We say that a link £ v is feasible in S if Eqn. Q] 
is satisfied for £ v . A set S is feasible if each of its link is feasible. 

Note that what we described above is the abstract SINR model. In the more commonly studied 
geometric SINR model, G vw is a polynomial function of d(s w ,r v ), where d(x,y) is the distance 
between two points x and y. Our results naturally apply to that model as well. Given that the 
geometric SINR model does not capture obstacles, reflections and other real life distortions, it is 
interesting to see what can be proven in the abstract model. 

Our setting where the powers are given as part of the input is often called the fixed power case, 
as opposed to the power control case where the algorithm can choose the power assignment. So 
far, research on fixed power has focused on oblivious power assignments, where the power of a link 
is a (usually simple) function of the length of the link [131 El USB 02] • Recently, a constant factor 
approximation algorithm to find the capacity in the power control case has also been achieved [19] • 
Unfortunately, none of these techniques appear to extend to the case of arbitrary fixed powers (for 
either arbitrary or geometric gain matrices). Yet, the problem of arbitrary fixed powers is not 



only natural, but has practical relevance, as commercial hardware often do not have the capacity 
of choosing precise powers to implement either an arbitrary assignment a la [19J, or to implement 
many of the oblivious power assignments found in literature. 
In this paper, we prove the following theorem. 

Theorem 1. Assume (L,P,G) is an instance of the capacity problem in the abstract SINR model, 
such that \OPT\ > i(l + e)\L\ for some e > ; where OPT is the maximum feasible subset of L 
using P. Then there is a polynomial time randomized algorithm to find a feasible set of size 0(e|-L|), 
with probability 1 — o(l). 

We do this by means of a semi-definite programming relaxation, which we show how to suc- 
cessfully round if the condition \OPT\ > ^(1 + e)\L\ holds. In addition, we discuss numerical 
experiments we have performed. These experiments show that the algorithm appears to work quite 
well on random instances, even better than the guarantees of Thm. [TJ 

Semi-definite programming has been a staple in designing approximation algorithms for NP-hard 
problems ever since the seminal work of Goemans and Williams on the Max-CUT problem [7]. It 
is interesting to note that the discrete "classical" problems closest to wireless capacity, namely 
the independent set problem and the graph coloring problem, have been fruitfully studied using 
semi-definite programming [15} 118]. The vertex cover problem, also relevant via its connection to 
the independent set problem, also has SDP-based approximation algorithms 117], Given this 
background, one may expect some of the techniques to easily carry over to the capacity problem. 
Yet that does not appear to be the case, at least not in a straightforward manner. A study of 
the aforementioned papers reveal that the discreteness of the problem plays an important role in 
the bounds. For example, in [18], the analysis proceeds by bounding the probability of vectors 
representing edges not being cut by a random hyperplane. Given the additive nature of the SINR 
model, it is not obvious how to extend that analysis to this case. There have also been a number 
of results for these problems on hypergraphs [21j [21 [3]. Though hypergraphs appear to be closer 
in spirit to the additive wireless model, they are still different, because the effect of each node on 
any other node doesn't change in the SINR model (as opposed to in a hypergraph, where it can be 
different based on which edge they are in). Thus, the (sophisticated) methods on hypergraphs do 
not appear to translate immediately to the SINR model either. Our SDP relaxation and rounding 
algorithms are quite simple in contrast to some of the previously mentioned work. Whether or not 
advanced techniques can be extended to the SINR model remains to be seen. 

1.1. Related Work. Moscibroda and Wattenhofer [24J were the first to study of the scheduling 
complexity of arbitrary set of wireless links. Early work on approximation algorithms produced 
approximation factors that grew with structural properties of the network [27 } I25| . IT] . 

The first constant factor approximation algorithm was obtained for capacity problem for uniform 
power in |8j (see also [T3]) in R 2 with a > 2. Fanghanel, Kesselheim and Vocking [6] gave an 
algorithm that uses at most 0(OPT + log 2 n) slots for the scheduling problem with linear power 
assignment P v = d(s v ,r v ) a , that holds in general distance metrics. 

Kesselheim obtained a 0(l)-approximation algorithm for the capacity problem with power con- 
trol for doubling metrics [19J. Around the same time, the first constant factor algorithm for all 
sub-linear, length monotone power assignments was achieved on general metrics [12j . Other re- 
cent studies in the SINR model include work on topological maps |16j . distributed algorithms for 
scheduling [11] . distributed power control [I] and auction based spectrum allocation |23j. 

2. SDP-BASED ALGORITHM. 

First, some notation. Vectors are denoted by x, etc. The standard 2-norm of the vector x is 
||x||. The i th entry of x is x(i). The inner product of vectors x and y is denoted (x ■ y). Define 
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9w = P v Gvv — j3N and g vw = P W G VW for v ^ w. Note that we can assume without loss of generality 
that g vv > 0, Vt>. Let OPT be a feasible subset of L of maximum size. Note that n = \L\. 
Consider the following program. 



max 2 (st ■ s), subject to 

V 

(s v ■ s)g vv >P\^2(s 

{s v ■ S) > 0, Vf 
(s„ ■ s^,) > 0, Vv, u; 

(4 ■ Sw) > («t> • + (sw • s) - 1, Vu, w 
|js^,|| 2 = 1, Vt> and ||s|| 2 = 1 . 

where s* v , s € M ra+1 . Each link has a vector variable s~l associated with it. The dot product of 
with a vector s denotes the (fractional) extent to which £ v is selected in the solution. 

Since the objective function and constraints are all linear functions of vector inner products, this 
problem is a SDP. Thus the program can be solved up to an additive error of e > in time that 
is polynomial in n and loge |28| . Since e can be made small enough to not matter, we will simply 
assume that the problem can be solved exactly. 

We can rotate the vectors to fix s = {1, . . . 0}, thus the above program is equivalent to: 



max s ^>(l)i s -t- 

V 

(2) s-* v (l)g m >/3\J2( 

(3) Se(l)>0,V« 

(4) {s v ■ sZj) > 0,Vv,w 

(5) (s* v ■ Sw) > + - l,Vv,w 

(6) ll s "i>l| 2 = l,Vu . 

Let us verify that this program is a relaxation of the maximum capacity problem. 
Lemma 2. The SDP is a relaxation of the original problem. 

Proof. Consider any optimal solution OPT to the capacity problem. For all £ v € OPT, set s~l 
s = {1,0,0,0. .. 0}. If £ v G L\OPT set 



1 if i = v + 1 
otherwise 



In other words, we make sure that each unselected link chooses a different position for the single 1 
in the vector. 

Given these assignments, Equations El H] and [6] can easily seen to hold. 

To show that Eqn. [2]is satisfied, first assume £ v G OPT . The following observation is immediate; 

Claim 2.1. Ifl v ,£ w G OPT then s^l) = s^(l) = (s~ ■ s„) = 1. Ifl v eL\ OPT then s* v (l) = 
and (s~l • sZj) = for any £ w ^ £ v . 
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Since l v G OPT, 

&v\±)9vv — 9v 



And, 



P I ^ {sy ' Siv)9vw 



v€OPT\{v} J \w<EL\(OPTU{v}) 



= P I X] 5ot 

\t«eOPT\{i)} 

where the second equality follows from the claim above. 

Now, since £ v G OPT, g vv > f3 (y2 we oPT\{v} 9wi\ (by Eqn [I]). Thus, the above two equations 
show that Eqn. [2] is satisfied when £ v G OPT . The case where £ v G" OPT is similar. 

For Eqn. [5l the following observations suffice: 

• If £ v ,£ w G OPT, fa . S %,) = 1 = 8- v (l) + s%,(l) - 1 

• If £ v ,£ w G" OPT , they have Is in different positions and (s^ • s^) = > + — 1 

• If £ v G OPT, £ w G" OPT, they have Is in different positions and (s v ■ s w ) = = 1 + 0-1 

□ □ 

Now we present our algorithm and the proof of Thm. [TJ We need two related definitions. Let 
8 V = max{s7>(l) — i, 0} for all £ v G L. Further, define L + = {£ v G L : S v > 0}. The algorithm is as 
follows. 

Algorithm 1 Capacityl 



Solve the SDP 

Select each link l v G L + with probability in to a set R 
Output {£ v G R : £ v is feasible in R} 



Lemma 3. If \OPT\ > (1 + e)n/2, £/ien E^eL+ ^« — T' 

Proof. Since |OPT| > (l+e)n/2, it follows that s«(l) > (l+e)n/2 (since the SDP is a relaxation 
of the original problem). Now by definition of 5 V , 6 V 

+ \> s v {l). Thus, 
E (^ + >(l + e)n/2 
=> 53 > (1 + e)n/2 - [L|/2 = (1 + e)n/2 - n/2 = y 

Observing that 5 V = for G" L + completes the proof. □ □ 

We can now prove the main Theorem. 
Proof, of Thm. [T] 

Assume that the random binary variable X v describes whether or not t v G L + is chosen into R. 
We observe that E(X„) = k^, according to the algorithm. 



Then for any £ v , 

E [ P [ J] 5«™ J J = E | /? | Y, 9v W X w 

K wER\{v} J J \ \w£L+\{v} 

p\ y 9v W nx w )\=pi y 9vwy 

(7) E [ p [ Y 9v W I ) = \p [ Y 9v W S v 



Now, by Eqn. El 



s v {l)g vv >P\ Y( s )9vw , Vi> G L + 



Since 5^(1) ^ ^ for v G L~^~ and (s^ • s w *)c/ vw is always non- negative, we get for £ v G L + , 

\w£L+\{v} 

\wEL+\{v} J \ w( zL+\{v} 

(8) > P J] <W™ 

\«>eL+\{i>} 

where the second inequality follows from Eqn. El and the first equality follows from observing that 
S v = s* v (l) - \ for £ v G L + . 
Then, for £ v G L + , 

F(£ v is infeasible in R) = P I /? 5^ J > <7 OT 

, n . ^ e (^(x:^ 6jR \ { „ } 1 

l y J s s r 

The first equality is the definition of infeasiblity. The first inequality is Markov's inequality. The 
last inequality follows from Equations [7] and El 
Now the expected size of the output is 

E G R : £ v is feasible in R}\) = Y p (4> G R and is feasible in R) 

e v eL+ 

= Y p (4> is feasible in J R)P(4 G R) 

tv€L+ 

-Z^22 4 ^ " - 8 
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21 26 31 36 41 

Max feasible set 

Figure 1 . OPT vs the average size of the set found by the SDP algorithm. In each 
case n = 61. 

The second equality follows from the independence of the events concerned. The first inequality 
follows from Eqn. [TUl The last inequality follows from Lemma El Thus the expected size of the 
feasible output is O(en). It is not difficult to boost the probability of getting a Q(en) size subset 
to complete the proof of the theorem. □ □ 

3. Numerical Experiments 

We ran simulations to test how well the algorithm does in practice. We used CVX, a package for 
specifying and solving convex programs using MATLAB [9]. We ran it on version 7.8 of MATLAB 
running on a Macbook with a 2 GHz Intel Core 2 Duo Processor and 2 GB of RAM. 

We generated a number of problem instances where n = 61 and \OPT\ = 21,26,31,36 and 41. 
The instances were generated as follows. To generate the feasible subset a large random instance M 
of links on the 2d plane was generated. Each sender s v = (s v (x), s v (y)) is a random point in a 450 x 
450 box. The receiver r v is defined by (s v (x) + random,, (2;), s v (y) + randomly)) where random^(x) 
and randomly) are sampled uniformly at random from [—20,20]. We generated corresponding 
gain matrices using the geometric SINR model setting a = 2.5 (thus G vw = Ti — - — We used 

both uniform (P v is a constant) and mecin power assignments — \\$v — 

II ) t° generate the 

gain matrix. We set the noise N = throughout the experiments. 

To generate the input instance G (which is a n x n matrix), we combined a subset of M with 
random entries. More specifically, first we retrieved a random feasible subset R of M (found 
greedily). This defined a R x R submatrix of G. The remaining entries were chosen iid randomly 
from [0, ft], where k was chosen large enough so that the remaining n—\R\ links would not contain 
a large feasible subset, thus R would be OPT for the instance. 

Though computationally slow (for n = 60 the SDP took a few minutes to be solved), the 
algorithm performed extremely well. Indeed, it took some time to come up with instances where 
the algorithm didn't have a perfectly integral solution. If the random entries of G corresponding 
to L \ OPT were too large (corresponding to a large k, meaning that L \ OPT contained only very 
small subsets that were feasible) or if OPT was too loosely feasible (ie, Eqn. [2] was far from being 
tight for most of the links), the algorithm did exceedingly well. 
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21/20/20 26/20/15 31/30 36/25 41/20 51/10 

Configuration 

Figure 2. OPT vs the average size of the set found by the SDP algorithm. In each 
case n = 61 and the links in OPT were generated using mean power. The labels in 
the x-axis describe the configuration of the instances. Thus, in the first case, the 
instance is an union of 3 feasible sets of size 21, 20 and 20, respectively, where the 
latter two are copies of subsets of the first one. 

Even after trying to make the problem more difficult, the algorithm did quite well, only degrading 
when OPT < n/2, for which we claim no theoretical guarantee anyway, though even in these cases 
the output was not unsatisfactory. Indeed, in all these cases, using the simple filtering (s* v ■ s) > 0.51 
identified OPT almost exactly. Our sampling algorithm, by design cannot achieve better than a 
factor 2 approximation in general, and that is almost what we achieved in all cases, as illustrated 
in Figure Q] for uniform power (the results for mean power were essentially identical) . 

As we have mentioned, in the above experiments, the algorithm sharply identified OPT. To 
create more ambiguous instances, we also tried the following. In this setting, we took a feasible set, 
and added copies of subsets of it. Thus the instance would be of the form L\ U L2 or L\ U L2 U L3 
where L\ is feasible, and L2,L% are copies of subsets of L\. One expects the solution to be more 
"spread out" in this case, and that is exactly what we found. The algorithm still performed rather 
well, even below theoretically guaranteed levels, though the behavior is somewhat different. Figure 
[2] demonstrates the case for mean power. 

4. Conclusion 

We have shown how to use semi-definite programming to approximate the wireless capacity 
problem in cases where the capacity is known to be large. It is an interesting question whether or 
not these results can be further improved, potentially using the power of geometric SINR model. 
Questions about the integrality gap and hardness of the problem (apart from what is known via 
the fact that the problem generalizes max independent set) also deserve attention. Though we have 
performed some preliminary numerical experiments, the efficacy of this method both in terms of 
accuracy and computational efficiency also is an interesting avenue of further investigation. 
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